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ABSTRACT 

The X-ray emission from X-ray binaries may originate in flares occurring when mag- 
netic loops anchored in the disc reconnect. In analogy with our Sun, Ha emission 
should arise as the accelerated electrons thermalize in the optically emitting disc, per- 
haps leading to correlated variability between X-rays, Ha and the optical continuum. 
We present simultaneous X-ray and optical high speed photometry of the neutron star 
low-mass X-ray binary Cyg X-2 to search for such correlations. The highest time reso- 
lution achieved is 5 ms in white light and 100 ms with a 3 nm filter centred on Ha. We 
find power on timescales > 100s (fiickering) in optical with a total r.m.s. of a few %, 
about an order of magnitude less than that seen in X-rays. We do not find significant 
correlations between the X-ray and optical fluxes on short timescales, hence cannot 
conclude whether magnetic flares contribute significantly to the optical emission. 

Key words: accretion, accretion discs — binaries: close — stars: individual {V1341 
Cyg) — X-rays: individual (Cyg X-2) 



1 INTRODUCTION 

■ Observations of low mass X-ray binaries (LMXBs) show 
\ that strong emission at X-ray energi es < 100 k eV is 
I nearly ubiquitous in these systems (e.g. |Do3|20o3)- The 
standard optically thick, g eometrically thin accretion disc 
iShakura &: Sunvaevlll973) cannot explain X-ray emission 
beyond a few keV. Thus, the observed hard X-rays as 
well as strong evide nce for extended "coronal" emis sion 
in eclipsing LMXBs (Ivan Paradiis fc McClintock|[l995l) re- 
quire more elaborate accretion flow models. A popular 
scenario, inspired by the solar corona, suggests that the 
cool Shakura & Sunyaev disc is sandwiched between two 
layers of hot plasma heated by reconnection of magnetic 
field l oops anchored in th e disc llGaleev. Rosner fc Vaianal 
Il979l: iHaardt fc Maraschi Il99ll) . The presence of such 
a magnetically-heated corona would not be surprising, 
considering that angular momentum in discs is most 
likely transported outwa rds by turbulent magnetic stresses 
iBalbus fc Hawlevlll99^ . In the context of this model the 
observed X-ray spectra are interpreted as time-averaged 
emission from multiple coronal flares, while the character- 
istic X-ray time variability seen both in neutron star and 
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black hole systems (e.g. . Ivan der Klij|2OO0l) is then related 
to a distribution of flare strengths and durations. 

In our Sun, bursts of hard X-ray emission from coronal 
flares are accompanied by simultaneous Ho? line emission 
when the beam of electrons accelerated during recon nection 
is thermalised in the photosphere (e.g.. IDu1M|1985() . If the 
origin of X-ray emission is similar to that in solar flares, cor- 
relations between Ha and X-ray emission should occur on 
timescales ranging from the rise time of the flares (~ O.Ols) 
to their life time (either the dynamical time of the accre- 
tion disc at the radius where a flare occurs, or the thermal 
time needed to dissipate the energy of the wound-up mag- 
netic field, up to ~ tens of seconds at the outer edge of the 
disc) . In the optical spectra of LMXBs, Ha is the strongest 
line (EW ~ 1-10 A, FWHM - 1000-3000 km s"^) and its 
double-peaked emiss ion profile is an unambiguous sign ature 
of the accretion disc Jvan Paradiis fc McC^ntockll995^ . The 
disc is optically thick, so the line has to be produced in a 
thermally inverted layer, presumably located above the pho- 
tosphere. However, if we pursue the analogy with the Sun, 
the heating source of this disc chromosphere, be it coronal 
flares, viscosity, reprocessing of X-rays or MHD waves, has 
yet to be understood. 

A significant fraction of the optical flux in persis- 
tent LMXBs comes from the reprocessing in the disc 
of high energy photons produced in the inner parts 
of the accretion disk, or in the NS boundary layer 
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Jvan Paradiis fc McClintoc^ Il995l: lO'Brien et all I2OO2I) . 

This effect is most clearly demonstrated by the optical 
echoes of NS X-ray bursts observed in some systems. A 
few second delay between these optical echos and the X-ray 
bursts is inter preted as the light tra vel time to the reprocess- 
ing site {e.g., lO'Brien et al.ll2002r . This delayed heating of 
the disk by non-local radiation can in principle be distin- 
guished from quasi-instantaneous local heating by particles 
accelerated in flares using high time resolution simultaneous 
X-ray and optical observations. 

While the X-ray timing properties of LMXBs are well 
established, little is known about short-timescale non-orbital 
variations in the optical. On timescales > 1 min flickering in 
white light (characterised by erratic flux changes by factors 
< 2) is a common occurrence in both persistent and tran- 
sient LMXBs. The NS LMXB Sco X-1 shows intermittent 
correlat ed X-ray /optical flar es related to non- local repro- 
cessing (|llovaiskv et alJ llQSCll: IPetro et aLlllQsill. LMC X - 
2 also displays similar behaviour jMcGowan et alJ l20o3) . 
Recent monitoring of low-luminosity black hole LMXBs 
has revealed intriguing optical flares on timescales of tens 
of seconds (Hvncs et al. 2003b). On the other hand, on 
timescales < 1 min, the variability of the optical contin- 
uum or Ha emission in LMXBs and its correlation with 
X-rays are still largely unexplored; such studies are ham- 
pered by the faintness of most systems, even when in out- 
burst. Only two LMXBs have well established X-ray /optical 
continuum correlated behaviour o n very short times cales 
(down to a few msl: GX 339-44 llMotch et alJ 1198,1) an d 
XTE J1118-F480 ijKanbach et al.ll200ll : iHvnes et al.ll2003a^ . 
In both the variable optical flux has been interpreted as syn- 
chrotron emission from energetic flares (^Fabian et aliil982t 
Idi Matteo. Celotti Fa.hia.nlll 997^. 

In this paper we report on simultaneous X-ray and 
Ha/white light high time resolution observations of Cyg X-2, 
an optically bright (1/ ~ 14.8) persiste nt LMXB with a neu- 
tron star primary (e.g.. lOrosz fc Kuulkers,. 1999^1 . Cyg X-2 
has rapid (hour-days) X-ray spectral variations closely re- 
lated to the timing properties: the slope/amplitude of the 
power law and the frequency of quasi-periodic oscillations 
(QPOs, f ~l-20 Hz) se en in the power spectr um vary with 
the X-ray colours (e.g.. Iwiinands et al.lll99^ . The optical 
flux varies by < 1 mag on timescale of hours to days, and 
flickering is reported in the brightest state. An ellipsoidal 
modulation on the orbital period of 9.8 days is observed with 
a V band amplitude of 0.1 mag. The donor star contributes 
about 50-70% of the total optical flux. Ha is known to vary 
on timescales of hours to days (see §4), but no systematic 
study of the Ha line variability has ever been undertaken. 
Here, we describe our search for correlated rapid variations 
between the X-ray and optical fluxes of Cyg X-2. 



2 OBSERVATIONS AND DATA REDUCTION 

We studied the variability of the optical continuum. Ha line 
emission and X-ray flux of Cyg X-2 at frequencies < 100 Hz 
using simultaneous data obtained on 2-4 August 2001 (UT) 
at the Hale 200-inch at Palomar and with the Rossi X-ray 
Timing Explorer {RXTE). 



2.1 X-ray data 

Cyg X-2 was observed with RXTE/PCA llSwank et alll99d) 
on three consecutive nights (2-4 Aug 2001), for three con- 
secutive orbits on each night. Data were analysed using 
FTOOLS v5.0. For all nine separate observations, we used 
data types collected by the Experiment Data System (EDS) 
in three modes (SB_125us_0_13_ls , SB_125us_14_17_ls , 
SB_125us_18_23_ls), which collected counts of 2.0-5.7 keV 
(centroid energies), 6.1-7.4 keV, and 7.8-9.8 keV. We made 
use of event mode data (E_125_64M_24_ls) in the 9.5-20 keV 
energy range. The time-tags associated with the counts were 
corrected to the solar-system barycenter in the DE405 time 
system using axbary. 

We extracted a spectrum from the Standard 2 EDS 
mode data from the first night's observation. Examination 
of this spectrum showed that the countrate in excess of the 
background by a factor of 2-3 orders of magnitude up to 
10 keV, and by a factor of 2-10 in the 10-20 keV energy 
range. We neglected the background for the spectral anal- 
yses herein. The 2.0-9.8 keV lightcurves (rebinned on 32 s 
from the original resolution of l/256s) for each of the 9 visits 
are shown in the bottom panel of Fig. 1. The dip in coun- 
trate during the third visit is typical of the flaring branch 
(see §3.1). 



2.2 Optical data 

High-speed optical photometry of Cyg X-2 was obtained us- 
ing a CCD in a multiple frame-transfer observing mode. The 
detector used is a SITe SI-502A thinned, AR-coated, back- 
illuminated 512x512-pixel CCD. The F/9 Cassegrain focus 
of the Palomar 200- inch telescope illuminated a 200x4 arc- 
sec slit, which was re-imaged onto the uppermost rows of 
the CCD. The image was demagnified to a plate scale of 
0.4 arcsec/pixel, yielding a slit size of 500x10 pixels. The 
un-illuminated portion of the CCD was used as storage, al- 
lowing 50 images to be accumulated before reading out the 
entire CCD. The transfer time, the time required to shift an 
image from the illuminated region of the CCD into storage, 
is short, so the shutter remained open during the accumu- 
lation of all 50 images, and was only closed during readout. 

The timing of the images in a single frame was con- 
trolled by a GPS receiver and time-code generator. The GPS 
timing allowed control of the start time (and duration) of 
each image to an accuracy of 1 /.is. The transfer time is 420 
/is (42 /is per row shift), which occurs immediately follow- 
ing the start of a given image. The effective exposure time 
is the same for every image on the CCD except for the first 
and last few images. The first image or images have an ef- 
fective exposure time that is reduced by the time it takes 
the shutter to open, which is approximately 13 ms. The last 
image has an effective exposure time that is increased by 
the time it takes the shutter to close, although the readout 
takes place during this time, resulting in some smearing. For 
this reason, the first image or images and the last image are 
discarded in the photometric analysis. 

CCD pixels can be digitized in 22.4 /is, or discarded in 
2.4 /is. To minimize the time spent reading the CCD, two 
windows are defined, which mark ranges of columns in which 
the pixels are digitized, with pixels in all other columns dis- 
carded. Each window was 40 pixels wide, which allowed the 
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Table 1. Log of Cyg X-2 Optical Observations. 



Date 


Band 


At 


start 


stop 


Aper 


S/N 


samples 


# frames 


samples 


2001 




ms 


hrs 


pixels 




per frame 




kept 


Aug, 2 


wl 


5 


4.107 


4.943 


7x7 


17.2 


47 


1200 


56237 


Aug, 2 


Ha 


100 


5.842 


8.173 


3x3 


5.3 


48 


1100 


52573 


Aug, 3 


Ha 


100 


29.490 


30.325 


3x3 


4.3 


48 


400 


19164 


Aug, 3 


wl 


10 


30.970 


34.801 


7x7 


10.5 


47 


4400 


113572 


Aug, 4 


wl 


5 


53.211 


54.193 


7x7 


16.5 


46 


1200 


50484 


Aug, 4 


Ha 


100 


54.759 


58.980 


3x3 


4.6 


48 


2000 


63041 



Ha refers to a 3 nm filter centred on 656.2 nm, wl refers to observations taken in white light; At is 
the time resolution of the data; start and stop times are in hours from MJD 52123.0 UT ; Aper. is 
the size of the (square) aperture used for photometry in pixels; S/N is the observed signal-to-noise 
of each At exposure; samples per frame is the number of exposures on each CCD frame; last column 
is the number of flux measurements kept after bad data points were removed (see §2.2). 
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Figure 1. Simultaneous optical (top, Palomar 200- inch) and X-ray (bottom, RXTE) lightcurves of Cyg X-2. The time is in hours from 
MJD 52123. Each point is a 32 s (elapsed time) average. The type of photometry (white light or Ho filter) is indicated above each optical 
lightcurve. Amag is the differential magnitude between Cyg X-2 and the comparison star. The X-ray count rate was derived from PCA 
channels 0—23 (2.0-9.8 keV). Cyg X-2 was on the normal branch during the first two X-ray visits, on the fiaring branch on the third and 
on the horizontal branch in the subsequent observations. Typical error bars for the 32 s photometry bins are 0.002-0.003 mag in white 
light, 0.014-0.018 mag in Ha and 1.2-1.5-10"^ log(counts/s) in X-rays. 



entire CCD to be read out in 1.5 s. Overhead in the com- 
puter operating system added up to one second to the total 
time required for each frame. The exposure times used in 
these observations were 5 ms, 10 ms, and 100 ms per image, 
giving duty cycles of approximately 10%, 20%, and 65%, 
respectively. 

The long slit, with two windows defined, ensured 
that each image included a bright comparison star. The 
comparison star used is t he photometric standard 35 of 
iHenden fc HonevcuttllT997l {V = 13.498 ± 0.001, B - V = 
1.194 ± 0.004) situated ~ 1' NNW of Cyg X-2. This star 
has been extensively monitored, and has shown no intrinsic 
variability down to 0.001 mag. Roughly half of the data was 
taken with a 3 nm FWHM filter centred on Ha and the re- 
mainder in white light (see Tab. 1). The seeing over three 



nights ranged from 0.8 to 1.2 arcsec FWHM in Hq, and from 
1.2 to 1.6 arcsec FWHM in white light. 

Lightcurves were constructed for each observation se- 
quence using square photometric apertures ranging in size 
from 1x1 to 11x11 pixels (0.4x0.4 to 4.4x4.4 arcsec). While 
the use of square photometric apertures, rather than circu- 
lar apertures, is not conventional, there is little difference 
for apertures composed of a small number of pixels. For the 
larger apertures, there was no measurable difi'erence in S/N 
when circular apertures were used, so for consistency, all 
photometric apertures used were square. 

For each image the aperture is placed at the centroid of 
Cyg X-2, known to a fraction of a pixel. The comparison star 
aperture is placed at the (fixed) mean distance from Cyg X-2 
to compensate for aperture losses. We checked that the PSFs 
of the two stars were similar and not affected by instrumen- 
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tal effects. Because every image includes botfi Cyg X-2 and 
the nearby comparison star, observed simultaneously under 
identical observing conditions, differential photometry al- 
lows the removal of temporal photometric throughput varia- 
tions. We do not correct for extinction hence the broad band 
white light data can be subject to an airmass-correlated vari- 
ation due to the colour difference between Cyg X-2 and the 
comparison star. However, since all our observations were 
done at airmass <1.2, with two-thirds of the data taken at 
airmass < 1.1, we feel this will have a negligible influence 
on the short timescale variability study. 

We ignored images for which the centroiding algorithm 
had failed, a cosmic ray had been tagged inside the aperture, 
or (more commonly) those images in which the photometric 
aperture extended outside the unvignetted slit edges. The 
background in each window of each image was calculated as 
an average over the window, for pixels more than 10 pixels 
(4 arcsec) from the window center. We then computed the 
signal-to-noise ratio of the flux from Cyg X-2 as a function 
of aperture. The maximum S/N was obtained for an aper- 
ture of 7x7 pixels in white light and 3x3 pixels with the 
Ha filter. The comparison star data was extracted using the 
same aperture and the flux ratio computed to correct for 
lost flux and atmospheric variations. The final lightcurve is 
therefore the differential magnitude between Cyg X-2 and 
the comparison star. The S/N, chosen aperture, total num- 
ber of measurements and number of measurements kept in 
each observation sequence can be found in Tab. 1. 

The lightcurves, rebinned with 32 s resolution, are 
shown in Fig. 1. Note that 32 s is an elapsed time and 
that the actual exposure time and S/N in each bin varies 
with the observational setup. The la error for each bin is 
obtained from the standard deviation of the averaged data 
points, hence implicitly assuming that the measurements are 
independent. Typical error values are given in column 4 of 
Tab. 2. 



3 ANALYSIS AND RESULTS 
3.1 X-ray data 

To characterise the overall behaviour of Cyg X-2 in X-rays, 
a colour-intensity diagram was plotted using 32 s flux aver- 
ages (Fig. 2). This diagram suggests that Cyg X-2 was on 
the normal branch during visits 1 and 2 (Aug 2), on the 
flaring branch during visit 3 (Aug 2) and on the horizontal 
branch during the visits of Aug 3 and 4. Power density spec- 
tra (PDS, see Ivan der K lis 1989) were then constructed for 
each visit using all the detected photons. The high frequency 
(1 Hz-128 Hz) part of the PDS is the average PDS of 1 s long 
segments with a binsize of 1/256 s. The low frequency part 
10~*-0.5 Hz) is computed from a 1 s rebinned lightcurve. 
The PDS are normalised to units of r.m.s.^ Hz~^, i.e., frac- 
tional variance per unit frequency, and the expected Poisson 
measurement noise bias was subtracted. No corrections were 
made for background and PCA dead time, which are negli- 
gible. 

The resulting PDS conflrm that Cyg X-2 moved through 
the different branches of the Z track during the observations. 
On Aug 2, the absence of detectable broad band variabil- 
ity at frequencies above 1 Hz except for a weak QPO at 
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Figure 2. Hardness-intensity diagram of Cyg X-2 on Aug. 2, 3 
and 4. The hard colour (y-axis) is the log of the count rate ratio 
between 5.7-9.8 and 2.0-5.7 keV (PCA channels 14-23 and 0-13). 
The intensity (x-axis) is from 2.0 to 9.8 keV. Each point is a 32 s 
average. The numbers indicate the position of the points from 
each X-ray visit (as labelled in Fig. 1). 



6±1 Hz during the second visit (1.3±0.2% r.m.s.) is typical 
of the normal/flaring branch. The Aug 3-4 PDS are typi- 
cal of the horizontal branch with strong broad band vari- 
ability peaking around 10 Hz and two QPOs varying be- 
tween ~24-34 Hz (3.1±0.2% to 4.4±0.2% r.m.s.) and ~50- 
63 Hz (2.2±0.3% to 2.3±0.3% r.m.s.). In addition, low fre- 
quency variability is present in all the observations at fre- 
quencies below 0.1 Hz. These results are in line with previ- 
ous, more exten s ive X- r ay timing studies of this source (e.g. 
KViinands ct al. 1997; Kuulkcrs. Wii nands fc van der KlisI 
11999: Wiinands fc van der Klis 2001.). We show in Fig. 3 
some of the logarithmically rebinned X-ray PDS obtained 
on different nights. 



3.2 Optical data 

The 32 s binned optical lightcurves (Fig. 1) show varia- 
tions above the measurement noise on timescales of minutes 
to hours. Quantitatively, the binned white light (resp. Ha) 
lightcurves have standard deviations of up to 0.04 mag (resp. 
0.06) compared to the measurement noise of ~ 0.003 mag 
(resp. 0.017; see §2.2). The values for each observation are in 
Tab. 2. Plotting the optical fluxes against the simultaneous 
X-ray fluxes, we find no obvious correlations between the 
two; but we note that the average white light flux was high- 
est (by about 0.15 mag, second night), when the X-ray count 
rate was also highest (close to the vertex of the horizontal 
branch) . 

PDS were constructed for the optical data using a slight 
modification to the standard technique used for X-ray tim- 
ing. For the optical data, we found the white noise level em- 
pirically instead of using Poisson statistics. The data were 
(fast) Fourier transformed and normalised by multiplying 
the amplitudes by 2/Var(/), where Var(/) is total variance 
of the data /. This ensures the powers follow a xi distri- 
bution when the lightcurve is dominated by (measurement) 



Simultaneous optical and X-ray high speed photometry of Cyg X-2 5 



white noise {i.e., each point is independent and randomly 
taken from the initial flux distribution). The optical power 
spectra typically showed power at low frequencies and white 
noise at high frequencies. In these cases we divided by the 
average power at the frequencies where measurement noise 
dominates to recover the xi statistics. The Fourier ampli- 
tudes give the fractional variance per frequency bin and, 
after subtraction of the constant white noise component, 
we plot the power x frequency in (dimensionless) units of 
r.m.s.^ to highlight the timescales with significant variability 
and ease comparison with the X-ray PDS (Fig. 3). 

The low duty cycle of the data acquisition scheme in- 
troduces large amounts of spurious power in a Fourier trans- 
form covering several frames. To circumvent this problem, 
we separated the calculation of the power spectrum into two 
parts so that the transformed lightcurves had duty cycles 
close to 100%. At high frequencies, we computed a power 
spectrum of the images in each frame separately, then av- 
eraged the power spectra of all the frames. Since there are 
between 46 to 48 images (photometry points) in each frame 
with a time resolution of 5 ms (white light) to 100 ms (Ha) , 
the resulting frequency range is about 5-100 Hz (white light) 
or 0.2-4.8 Hz (Hq), depending upon the observational setup 
(see Tab. 2). At lower frequencies we took the average flux of 
each frame (the sum of all images in a frame) and computed 
the power spectrum of the resulting lightcurve. The result- 
ing frequency range is about 10"^-0.05 Hz (see Tab. 2). The 
combined power spectrum is free of power leakage at the ex- 
pense of a loss in frequency resolution and coverage. At high 
frequencies, the averaging of the individual spectra diminish 
the measured power from a highly coherent signal, because 
the contributions from each frame are not summed in-phase. 

Significant power is present in all of the datasets at 
very low frequencies (below 3 x 10""^ Hz in white light and 
10"'^ Hz in Ha) except during the short August 3 Hq obser- 
vation. The total r.m.s. of the flickering is estimated by inte- 
grating over the frequencies where power is detected above 
the 3(T single trial level. We obtain values of 1-3% r.m.s. (see 
Tab. 2) which are consistent with the standard deviation of 
the 32 s binned lightcurve. The optical variability is about 
an order of magnitude smaller than that seen in the X-ray 
lightcurves on the same frequency range. 

At higher frequencies, statistically significant power is 
detected in the white light observation of Aug. 3, when 
Cyg X-2 was on the horizontal branch with a high X-ray 
flux. The total r.m.s. is about 1% (0.02-0.07 Hz). Unfortu- 
nately, this is the only detection of optical variability at high 
frequency in our limited dataset. Independent confirmation 
would be desirable. The very low frequency power in both 
optical and X-rays is reasonably approximated hy Py 
{i.e., vPv ~ in Fig. 3). We obtained an upper limit to 
the flickering power in the cases where none was detected 
by calculating the total r.m.s. of power with this particular 
shape that would have yielded a 3(t detection (typically a 
few %, see Tab. 2). We also computed 3cr upper limits to the 
r.m.s. of a sinusoidal signal. Both of these upper limits take 
into account the number of trials in the frequency range. All 
of the results may be found in Tab. 2. 



3.3 Cross-correlation 

We cross-correlated each of the nine simultaneous X-ray and 
optical datasets (Tab. 1). The X-ray data, which have the 
highest time resolution, were corrected to the local Palomar 
time and rebinned in an identical way to the optical data set, 
yielding two continuous lightcurves with the same resolution 
(At = 5, 10 or 100ms). T he cross-correlation fun ction at a 
lag kAt is defined as (e.g. lEdelson fc Kroliklll98^ : 

c7(feAt) = i- ;^ (1) 

M ^-^ (TXCTopt 
j — i—k 

where o is the optical lightcurve, x the X-ray, the bar sym- 
bolises the average of the dataset and o are the standard 
deviations. C(fcAt) is an average over the M pairs of points 
which are fcAt apart. The variance of the M averaged dis- 
crete correlations gives an estimate of the error on C. The 
above is efficiently computed in Fourier space for continuous 
lightcurves with the same At: 

^(^^^^^^{ f'hf!:fu }.' 

where F is the (fast) Fourier transform, * is the conjugate 
and the w axe window functions (a timeseries with value 1 
when data are taken, elsewhere). 

There were no obvious trends on timescales >10 min 
where the datasets are not long enough to clearly estab- 
lish correlations. For the cross-correlations below 10 min, 
we subtracted piecewise linear flts in 600 s increments to the 
datasets to smooth out any underlying long timescale varia- 
tion. The X-ray and optical data did not show any signiflcant 
correlations down to 5 ms (100 ms). Cross-correlating with 
different energy selected X-ray lightcurves (2.0-5.7 keV, 5.7- 
9.8 keV, 9.8 keV and above) gave identical results. 



4 DISCUSSION 

Our dataset is consistent with previous studies reporting 
variations in the optical continuum of Cyg X-2 of ~0.5 mag 
in days, ~0.1 mag i n hours and ~0.05 mag in minutes 
teristiai^^lj 119671: iKruszewski et aP Il974 I iKilvachlovl 
19781: iBeskin et al.l [1972) • Orbital ellipsoidal modulation 
aside, there are few well-established features in the optical 
variability of Cyg X-2. There is a trend for the continuum 
to become bluer in U -B as Cyg X-2 becomes brighter {B- 
V staying ~ constant : lLvutvi fc Sunvaevlll973 : iBasko et alJ 
[1976). Varying fluxes and FWHM on timescales > 1 hr were 
also n oticed early on in Hq, H/3. Hell A4686 and th e Bowen 
blend djohnson fc Golsonlll969l : ICowlev et al.lll^79^ . 

The evolution from the horizontal to the flaring 
branches in Cyg X-2 like sources has been proposed to re- 
flect an in crectsing mass accre tion rate in the disc (see dis- 
cussion m iHoman et ai]l2002^ . This could be due to fluctu- 
ations in the mass transfer rate from the donor star or to 
fluctuations in the angular momentum transport processes, 
presumably on the disc viscous timescale. Interestingly, Hq 
is sometimes in absorption |Ch u vaey.l9_76-) ■ in dwarf novae, 
the Balmer lines evolve from emission to absorption during 
the rise to outburst when the mass a ccretion rate increases 
(e.g. ISzkodv. Piche fc Feinswo^liggoD . However, mass accre- 
tion rate changes are hard to reconcile with rapid (hours) 
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Table 2. Optical Variability of Cyg X-2 



Date 


Band 


At 


err* 


cr* 


^min ^max 


pulse 


total 


^min ^max 


pulse 


total 






ms 


10-3 


mag 


mHz 


% r 


m.s. 


Hz 


% r. 


m.s. 


Aug, 2 


Ha 


100 


14 


57 


0.12-66 


<0.4 


3.4 


0.21-5 


<1.2 


<2.5 


Aug, 3 


Ha 


100 


18 


24 


0.30-66 


<0.7 


<1.6 


0.21-5 


<1.9 


<3.8 


Aug, 4 


Ha 


100 


16 


37 


0.07-66 


<0.5 


2.1 


0.21-5 


<1.4 


<2.8 


Aug, 2 


wl 


5 


2.4 


14 


0.33-200 


<0.2 


1.0 


4.35-100 


<0.3 


<0.6 


Aug, 3 


wl 


10 


2.0 


41 


0.07-160 


<0.2 


2.2 


2.12-50 


<0.3 


1.1 


Aug, 4 


wl 


5 


2.6 


14 


0.28-160 


<0.2 


0.8 


4.35-100 


<0.4 


<0.7 



wl is white light; err* is the estimated error of a 32 s photometry bin (in mmag) and cr* is the standard 
deviation of the 32 s binned lightcurve (in mmag; see Fig. 1); the optical variability is studied in two frequency 
ranges (see §3.2 for details): t'min-'^max give the corresponding frequency ranges; pulse gives the r.m.s. upper 
limit for a sine signal; total is the broad band r.m.s. variability or upper limit (assuming F^, oc coloured 
noise). 



Aug. 2 [X-ray,Ha] 



Aug. 3 [X-ray, white light] 



Aug. 4 [X-ray, Ha] 




H 1 1 1 1 1 h- 



H 1 1 1 1 1- 



• . t 



10-3 10-2 10-1 1 10 loz io-> 10-= 10-2 10-1 1 10 loz 10-« 10-3 iQ-z lo-l i lo 102 
1/ [Hz] V [Hz] V [Hz] 



Figure 3. Temporal variability power spectra for the X-ray (top) and optical (bottom) data on different nights. White noise has been 
removed, the power normalised to units of r.m.s.'^ Hz— ^ and multiplied by frequency. The frequency range that contributes most to 
variability and its r.m.s. are thus easily identified. We also plot 3(t (single trial) error bars to help assess the statistical significance of 
the various features. The bottom of the first panel shows the power spectrum of the long Ha dataset taken on August 2; top shows the 
X-ray power spectrum from the simultaneous RXTE data (two visits, see Fig. 1). Second (resp. third) panel is for the long white light 
(resp. Ha) dataset taken on August 3 (resp. 4). The power spectra have been logarithmically rebinned. Downward arrows in the X-ray 
panels highlight the location of the QPOs discussed in §3.1. 



variations clianges unless they originate in the inner disc 
where tlie timescales are sliort. 

Variations in tlie inner region can afltect the opti- 
cal output of the disc via non-local irradiation. The ob- 
served correlation between liigli ionisation line fluxes (Hell, 
Clll/Nlll) and the optical continuum of Cy f^ X-2 does sug- 
gests that reprocessing plays a major role (Il gyaislcv et alJ 
Il979 : ^l,m Paradiis ct al. 1990; O'Brien 200 j). Some of 
the long timescale X-ray fluctuations of Cyg X-2 have 
been argued to be due t o varying obscuration, perhaps 
by an inner disc warp llKrmllte rs fc van der KlisI Jl99,4 
iKuullters. van der Klis fc VaugliaJl996l : IVrtilek et alJl99l 



|2Qfl^. Changes in the albedo or geometry of reprocessing 
would lead to complex correlations between X-rays and op- 
tical llEsin. Lasota. fc Hvne j l2000l) . O'Brien (2001) found 
that the optical and X-ray fluxes were correlated on the 
horizontal branch, but anticorrelated on long timescales 
in the normal and flaring branches. Curiously, on shorter 
timescales (~1 min), the X-ray and optical on the flar- 
ing branch were correlated. O'Brien proposed that several 
components with different intrinsic timescales contribute to 
the optical emission. Our observation of an order of magni- 
tude lower variability in optical compared to X-rays on long 
timescales (§3.2) is consistent with reprocessing of a fraction 
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of the X-rays, but we did not find significant correlations to 
support this. 

As argued in §1, variability and X-ray /optical correla- 
tions on short timescales could provide constraints on emis- 
sion sites in Cyg X-2. Previous studies found no variability in 
the optical continuum below ~ 1 min |Chevalier et'al]ll976l : 
lAuriemma et alJll978l : lllovaiskv et al.lll978D . With a factor 
2-5 better sensitivity and a longer dataset, we find some ev- 
idence for broad band variability at frequencies of few Hz in 
the continuum when Cyg X-2 was optically bright and on the 
horizontal branch; but we find no correlations with the X-ray 
noise component that peaks at ~ 10 Hz in th i s state . For 
rapid line variability, ICanizares fc McClintocl^ il975l) give 
an upper limit of 8% to the pulsed fraction of the Hen flux 
between 0.1 and 125 Hz. We do not flnd variability in the 
Hq line on timescales shorter than 1 min. 

In contrast with Cyg X-2, the other optically bright Z 
source, Sco X-1, shows correlated X-ray and optical flares 
on a 1-60 s timescale when it is on the flaring branch 
Jllovaiskv et al.lll98d : iPetro et al.lll98iri . The optical flares 
are 20 s filtered versions of the X-rays, probably as a result 
of non-local reprocessing. Similar behaviour is not seen in 
Cyg X-2, although the X-ray flux and variability are com- 
parable. A different irradiation geometry, as mentioned in 
the previous paragraph, might explain their absence. 

We did not find evidence for correlated X-ray - Ha vari- 
ations due to solar-like fiares. However, Ho? emission could 
be powered by fiares that are too weak to be detected indi- 
vidually and occurring at a fast enough rate that co rrelations 
are lost. For instance. IPoutanen fc FabiarJ il999ft discuss a 
model for the X-ray timing characteristics of Cyg X-1 in 
which the rapid variability is due to fiare avalanches in the 
disc. Only rare large fiares from the high luminosity end of 
the distribution function can be distinguished in the X-ray 
fiux of Cyg X-1 ICicrliriski & Zdziarski 200^. 

Calculating the thin disc respons e to irradiation by high 
energy electrons produced in fiares, IWilliams fc Maletestal 
i2002h flnd that an avalanche model can explain the H/3 
emission and the optical continuum flickering observed in 
cataclysmic variables (which have accretion discs compara- 
ble to those in LMXBs). Only 4% of the fiare energy is used 
to produce the variable V band fiux. A similar model might 
apply to Cyg X-2 with, perhaps, more of the fiare energy de- 
posited in the disc to take into account the higher optical/X- 
ray fiux ratio. The high number of active regions required 
is likely to have a significant impact on the X-ray spectral 
properties (e.g.. IdTMatteo. Celotti fc FabiarJ 119981) . X-ray 
emission by bremsstrahlung or Compton u pscattering is not 
modelled bv lWilliams fc Maletestal i2002f) . If flares give rise 
to the broad band X-ray noise (with r.m.s. levels of a few 
%) then our upper limit of ~ 3% r.m.s is too high to detect 
their Ha signature. 



5 CONCLUSION 

We obtained Ha or white light lightcurves of Cyg X-2 with 
high time resolution simultaneously with X-ray data. Two 
(mutually non exclusive) mechanisms may link the two wave 
bands: non-local reprocessing of X-ray photons emitted close 
to the compact object or local reprocessing of particles ac- 
celerated in magnetic fiares. We find broad band variabil- 



ity on timescales longer than a minute and some evidence 
on sub-second timescales in one white light observation. No 
correlations are found. Our observations cannot be used to 
distinguish between non-local or local reprocessing. The lack 
of correlations could be explained either by unfavourable ir- 
radiation circumstances (non-local case) or because large, 
identifiable magnetic flares are scarce (local case). 

The prospects offered by rapid optical photometry have 
proven difficult to fulfill. Part of the difficulty lays in as- 
sembling a large, coherent, simultaneous X-ray and opti- 
cal dataset with secure calibrations from which trends can 
emerge (i.e., the strategy behind the success of X-ray timing 
studies). More problematic are the small amplitude of the 
expected variability and limited number of bright enough 
objects for which sensitive enough observations and com- 
parisons can be made. In principle, the present generation 
of 8-10 m class telescopes could remedy that. 
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